Solution of the Multi-Channel Anderson Impurity Model 
— Ground State and Thermodynamics — 
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We present the solution of the SU(N)®SU(M) Anderson impurity model using the Bethe-Ansatz. 
We first explain what extensions to the formalism were required for the solution. Subsequently we 
determine the ground state and derive the thermodynamics over the full range of temperature and 
fields. We identify the different regimes of valence fluctuation at high temperatures, followed by 
moment formation or intrinsic mixed valence at intermediate temperatures and a low temperature 
non-Fermi liquid phase. Among other things we obtain the impurity entropy, charge valence, and 
specific heat over the full range of temperature. We show that the low-energy physics is governed 
by a line of fixed points. This describes non-Fermi-liquid behavior in the integral valence regime, 
associated with moment formation, as well as in the mixed valence regime where no moment forms. 



I. INTRODUCTION 

Heavy-fermion materials have been a source of interest 
and puzzlement since the experimental discovery of su- 
perconductivity in CeCu2Si2 in 1979. 1 By the turn of the 
1990s the non-Fermi liquid character of these materials 
was coming to the center of attention (for a recent review 
see Ref. 2] ) , coincidentally with the interest on marginal 
Fermi liquids generated by the normal state of high-T c 
superconductors. Since then, the number of examples 
of violation of Landau's Fermi liquid theory among lan- 
tanides and actinides has multiplied. 

Current theories trying to explain the non-Fermi liq- 
uid behavior in d- and f-electron metals can be classified 
into three broad categories: (i) models based on multi- 
channel Kondo physics, (ii) models considering the prox- 
imity of a quantum critical point, and (iii) models based 
on single-channel Kondo physics but in the presence of 
disorder that induces a distribution of impurity energy 
scales. These three ingredients are not mutually exclu- 
sive and a number of recent theories try to address their 
effects in different combinations (see |2( for references). 
In this article we will be concerned with the first class of 
models. This approach originated with the work of Cox^ 
in turn motivated by the unusual experimental results in 
the heavy-fermion compound UBC13. 4 He argued that the 
notably weak magnetic field dependence of the specific 
heat of this material excludes the usual magnetic Kondo 
effect and proposed instead that the observed anoma- 
lous behavior derives from the quenching of quadrupolar 
degrees of freedom. In this case, the spin of the con- 
duction band electrons plays the role of a channel degree 
of freedom. This is the two-channel quadrupolar Ander- 
son model, describing tetravalent uranium impurities in 
a cubic-symmetric metallic matrix. The model was later 
generalized to include other crystal symmetries as well as 
more complicated crystal field splittings; for a review see 



Ref. [f| • In more detail, Hund's rules and spin-orbit cou- 
pling in the presence of a cubic crystalline electric field 
lead to the modeling of a U ion in a Bei3 host by a 
Kramers doublet in a 5f 3 configuration and a quadrupo- 
lar (nonmagnetic) doublet T3 in the 5f 2 configuration. 
The doublets hybridize with conduction electrons in a 
representation carrying both spin and quadrupolar 
quantum numbers. This single-impurity approach was 
not uncontroversialAi For instance, whether the energy 
splitting between the two doublets is sufficiently large for 
a quadrupolar Kondo scenario to be viable is still largely 
unresolved. Aliev et al. presented experimental evidence 
suggesting that both doublets may in fact be nearly de- 
generate, pointing to a mixed-valent state with a novel 
type of interplay between magnetic and quadrupolar two- 
channel type screening. 8 

In this article we will study the general multi-channel 
Anderson impurity model that includes as particular 
cases the two-channel model and, to some extent, most 
of its generalizations alluded to above. 5 The so-called 
SU(N) <g) SU(M) Anderson impurity model, in its 
pseudo-particle representation, is given by the following 
Hamiltonian: 



#MchA = #host + £q ^ b^ba + E s ^ ftf<r+ 
a a 

+ VJ2{P*b^ a ,a (0) + <„ (0)&U 



subject to the constraint: J2 a b\b a + J2a faf" = 1- The 
first term in the Hamiltonian describes the host in which 
the impurity is embedded. For our purposes, we model 
it as a linearized Fermi band, 
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FIG. 1: Impurity level scheme for the case of an Anderson 
model with a three-fold degenerate magnetic configuration 
of energy e s and a five-fold degenerate quadrupolar one of 
energy e q . The energy difference e — e 
an arrow. 



e q is indicated with 



The second and third terms model two multiplets with 
energies e s and e q and quantum numbers a S SU (N) 
and a S SU(M), respectively^ We will refer to this 
two quantum numbers as generalized spin and flavor - in 
reference to the Kramers (magnetic) and non-Kramers 
(quadrupolar) doublets of the two-channel case. The last 
term in the Hamiltonian describes the hybridization of 
the host electrons with the impurity. 

An illustrative level scheme is given in Fig. ^ where 
the magnetic configuration is taken to lie lower in en- 
ergy that the quadrupolar one. The particular example 
is for an Anderson model with SU(N = 3)® SU(M = 5) 
symmetry. The impurity can switch between the two 
configurations by giving or taking an electron from the 
host band, a process that takes place with an overlap 
amplitude given by V. 

As a result of the constraint, the impurity Hilbert space 
is restricted to N + M states. Namely, N flavorless spin 
states \<j) = fl |0) plus M spinless flavor states \a) = 
b\ |0). Defining these states explicitly, we can rewrite the 
Hamiltonian in a different notation that automatically 
accounts for the Hilbert space restriction, 

-ffMchA = #host + £q ^ | a) (a\+e s J2\a) (a\ + 

a a 

+ v£[k> HVw (o) + VC H- 

a,er 

Both forms of the Hamiltonian are completely equiva- 
lent when the constraint is treated exactly and both are 
widely used in the literature^ It is the constraint acting 
in the Hilbert space that plays the role of a strong inter- 
action and renders the problem highly non-perturbative. 

As a side remark notice that if either M or N 
is put to one, the model reduces to a degenerate 
single-channel Anderson model in the limit of infinite 
Coulomb repulsion (a magnetic or a quadrupolar version 
of it, respectively) Ai£ While the standard single-channel 
SU(2) Anderson modelii was found to be integrableji2ii£ 
its SU(N) generalization is not integrable except for the 
strong repulsion limit when the impurity is constrained 
not to exceed single occupancy 14 i^ 



Even though the multi-channel Anderson model was 
put forward more than fifteen years ago, progress in its 
theoretical understanding has been slow. In the two- 
channel case (N = M = 2) most of the early knowl- 
edge of its unusual physics came from the integer va- 
lence limit. In this limit the model maps onto the two- 
channel Kondo model for which the Bethe-Ansatz solu- 
tion was available >iSiiLi& Also Numerical Renormaliza- 
tion Group (NRG) 1 ^^ and Boundary Conformal Field 
Theory (BCFT) 21 studies were carried out. Only more 
recently some progress was made in the study of the 
mixed valence regime of the two-channel Anderson model 
using Monte Carlo22 and NRG 23 i 24 methods. The more 
general multi-channel case is, however, not quite within 
the present reach of NRG and other approaches. On 
the other hand, the large N and M case constitutes 
the natural starting point for alternative approaches like 
the Non-Crossing Approximation (NCA)fSS* 2 ^ the con- 
serving slave boson theoryr^L or other types of 1/N- 
expansions. 28 Since the general multi-channel Anderson 
model can be regarded as the extension of the infinitely- 
repulsive degenerate single-channel Anderson model to 
the multi-channel case, the question about its integra- 
bility arises naturally. During the last couple of years 
the integrability was established, opening up the possi- 
bility of a full understanding of the model. In previ- 
ous work, we presented the Bethe-Ansatz solution for the 
two-channel case^ Subsequently, we developed the crit- 
ical low-energy theory of that model using BCFT and 
combining it with results from Thermodynamic Bcthc- 
Ansatz^i 



A. Preview of Main Results 

In the present work we will give a detailed and, to a 
large extent, self-contained account of the Bethe-Ansatz 
solution of the general multi-channel Anderson model. 
We shall show that, as in the case of the two channel 
model, the low energy physics of the impurity is governed 
by a line of boundary fixed points with a non-trivial resid- 
ual impurity entropy that is constant along the line. We 
shall identify two energy scales (Th and Tl) that gov- 
ern the screening process of the impurity degrees of free- 
dom. The screening occurs in two stages parameterized 
by these scales as will be seen, for instance, in the tem- 
perature dependence of the impurity entropy. The mixed 
valence regime will be discussed in detail, stressing not 
only the differences, but also the unexpected similarities 
with the integer valence cases. 

In Fig. |21 we present the picture that emerges for the 
different regimes of the model. As a function of tem- 
perature and energy difference between impurity config- 
urations (e = e s — e q ), we will characterize the different 
regimes (FV, LMM and LQM, and FP; see figure cap- 
tion). In particular, moment formation takes place (for a 
fixed and suitably large e) as the temperature falls below 
Tr; the moment being magnetic or quadrupolar depend- 
ing on the sign of e. It is then screened as T is further 
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FIG. 2: Schematic representation of the two temperature 
scales indicating the crossovers among different regimes: fluc- 
tuating valence (FV), fixed point (FP), and local magnetic 
and quadrupolar moment regimes (LMM and LQM). 



lowered below Xl when the system is governed by the 
infrared fixed point. Notice that the moment formation 
region becomes smaller as |e| is reduced and completely 
disappears at e ~ 0. In fact, over the whole mixed valence 
region (|e| < A = V 2 /2), valence fluctuations suffice to 
prevent moment formation. 

Another way of presenting the picture is to define two 
energy scales (T s and T q ) associated with the spin and 
quadrupolar degrees of freedom. These scales cross each 
other in the intermediate valence region and interchange 
roles as the high-temperature (Tjj) and low-temperature 
(II) scales that indicate the transition zones among dif- 
ferent regimes. The meeting of the two scales in the inter- 
mediate valence region signals the absence of a moment 
formation regime between the high and the low temper- 
atures: the fixed point is reached without prior moment 
formation. 

The rest of this article is organized as follows: in the 
next section we will discuss the scattering matrices and 
the integrability of the model (Sec. ITU; m the subsequent 
one we will present some necessary formal developments 
in the theory of quantum inverse scattering and the equa- 
tions derived from them (Sec. 1111(1 : in the following two 
we will discuss the thermodynamics of the model - we 
shall give formal derivations (Sec. llVfl and analytic and 
numerical results (Sec. [Vj -; in the last section we will 
provide a summary and a discussion of our main results 
as well as an outlook of their applications in the theory 
of heavy fermions fSec. IVI|) . 

II. S-MATRICES AND INTEGRABILITY 

The Hamiltonian -HmcIiA conserves the number of 
fermionic excitations (electrons), 

n c = J2 I (*) ^ (*) + E fU" 

a, a cr 



allowing us to study the system for an arbitrary but fixed 
value of N e . We take iV c = as our reference sector, an 
M-degenerate eigenstate with energy e q . Our strategy 
will be the usual one in coordinate Bcthc-Ansatz: we 
solve the system for N c — 1,2, .. . and then generalize 
the solution to arbitrary values of N e . Subsequently, as 
the N e — > oo limit is taken (discussed in later sections), 
the field theory is recovered. 

A. Electron-Impurity Scattering Matrix 

When there is only one electron present in the system 
(N c = 1), the most general one-fermion state has the 
following form: 

l fc ) = E / FL,p(x)€Ax)\f3)+G*\c-) . 

To determine the eigenstates of the Hamiltonian in this 
sector, satisfying H \ k) — E k \k), we apply the Hamilto- 
nian to this generic state and derive the first quantized 
Schrddinger equations in the sector. Using the expres- 
sion E k = k + e q for the eigenenergies of eigenstates |fc), 
we read off the equations: 

(-id x - k) F k a . p (x) + 6 (x) V 5?G k = 

(e-k)G k +^ a , <0 ,VS^F^. p , (0)=0 

where the energy difference e — e s — e q was introduced. 
We call F k a .p (x) and G k the wavefunctions. Eliminat- 
ing G k from the first equation and setting F* a .g (x) = 
e lkx F k a .p (x) we have (repeated indexes are summed 
over) 

(k e) {-id x ) F k aa .^ (x) + VH (x) (0) = . 

We make the following Bethe-type ansatz for the coor- 
dinate dependence of the F-wavef unction: F k a .p (x) = 

[(6(-x)I + 6(x)S)A]*, where S will be called the 
electron-impurity scattering matrix and A is an arbitrary 
vector in the internal space of the one-electron sector. In 
order to fully define the ansatz we adopt the following 
convention for the step function: 6(0) = 1/2. Let us 

introduce the operator [Q]".«f = 5^8^, that acts in fla- 
vor space and has the property Q 2 = M Q (recall M is 
the number of values that the index a takes) . Since A is 
arbitrary, one has a matrix equation for S. Its solution 
is 

Sl,0-Il,0- 2(fci _ e)+iW2 Ql;0 
e — iS(ki— e) _ ^ 

= Il ' + M Ql;0 ' 

where we used the index ' 1 ' for the only electron present 
in the system and introduced the use of the index '0' 
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for the impurity (the notation ko = £ will be also used 
latter). For the second way of writing S^o we introduced 

the phase 5 (k — e) — 2 arctan ^f^jy ■ 
It is easy to verify unitarity, 

SS^ = Si : oSo,i = I . 



B. Electron-Electron Scattering Matrix 

Consider now the case when there are two electrons 
present in the system (iV c = 2). The most general two- 
fermion state that we can write has the following form 
(all indices are summed over): 



\hk 2 ) = J J F*^ a>oa;/J ( Xl ,x 2 )^ aiai (xx)^ aaaa (x 2 ) \0) + J G k a \% !<72 ( Xl )^ aiai ( Xl ) \a 2 ) . 

Now again we apply the Hamiltonian in order to obtain the first quantized Schrodinger equations for eigenstates with 
eigenenergies E klk2 = k\ + k 2 + e q . We arrive at the following set of differential equations: 

E„=i, 2 (~id Xn - k n ) [AFtX a2 ^p ixuxt) + VA [S (x 2 ) 5^G k a \% a2 (a*)] = 
(-id Xl -k 1+ e~k 2 ) G k Q \% 0a2 ( Xl ) + V5(x 2 ) [AFfc£ aaaa . p (x u x 2 ) = 

I 



where A = I— V xsq = I— p^p^p? is (twice) the antisym- 
metrizer in coordinate, spin, and (quadrupolar)-flavor 
space, expressed in terms of P X ,P S ,P 9 the permutation 
operators that act in the spaces indicated. 

Now we make an ansatz for the F-wavefunction similar 
in spirit to the one we made in the one-electron sector. 
There are six regions in configuration space, correspond- 
ing to the six possible line arrangements of the impurity 
and the two electrons. The region with, say, electron 
'1' to the immediate left of the impurity is related via 
the S-matrix S^o to the region with electron '1' to the 
immediate right (the position of electron '2' remains un- 



changed). A similar role is played by S 2 ,o- However, we 
still need to determine how to relate the regions where 
electrons T' and '2' exchange places (away from the im- 
purity). 

Consider then, 

pfelfe2 ( \ _ ik 1 x 1 +ik 2 x 2 pfcifc 2 / T \ 

with the expression for F describing the six different ar- 
rangements of particles, 



^™;/3 fa. a*) = ^ M (~ x ^ 1 + M {-xx) S + 6 (-a*) 9 (x 2 ) S 2 , + 

+9 (-x 2 ) 9 ( Xl ) Si, S + 9 ( Xl ) 9 (x 12 ) S li0 S 2 , + 9 (x 2 ) (x 21 ) S 2 ,oSi l0 S) A 

I 



Here Si,o and S 2j o are the electron-impurity S-matrices 
found above and S = Si l2 is the electron-electron S- 
matrix that we seek to determine in this subsection. As 
above A is an arbitrary vector in the internal space of the 
two-electron sector, determining the state of two elec- 
trons and the impurity. To define unambiguously the 
ansatz we adopt the regularization 9 (0~) 9 (0 + ) = 0, 70 
consistent with the first order character of the differen- 
tial equations. Let us mention that this ansatz assumes 
that the same momenta k\ , k 2 characterize the wavefunc- 
tion in all six regions {i.e. orderings). This is at the heart 
of the ansatz and will be shown to be valid later when 
we discuss the Yang-Baxter conditions. 



Inserting the wavefunction into the first of the 
Schrodinger equations above we verify after some alge- 
bra that the equation holds, determining uniquely the 
form of the G-wavefunction. At this stage the electron- 
electron scattering matrix remains arbitrary and we turn 
our attention to the second Schrodinger equation. Car- 
rying out the algebra, we find that the equation holds 
provided the following matrix constraint on the electron- 
electron S-matrix is obeyed: 

(S 2|0 - I) (Si.oS - I) - P qs (Si, - I) (S 2 , - S) = 

where the matrices without indexes act on the internal 
space of the two electrons. A careful examination of this 
equation reveals the presence of an overall left-prcfactor 
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Q2 : o- Since this operator is not invertible, the solution 
of the constraint is not unique and there is still a certain 
amount of freedom left in the choice of S. 



1. An Integrable Solution 

The equation above does not have a unique solution for 
the two-electron S-matrix S. For our purpose we need to 
identify, however, a particular solution having the appro- 
priate physical properties that ensure the generalization 
and consistency of the ansatz to any number of electrons 
N c . These requirements are: 



Unitarity Condition 


Sij'Sji = I 


Locality of the Scattering 




Yang-Baxter Relation 





TABLE I: Three necessary conditions for integrability. 

The first of these relations assures the reversibility of 
the scattering processes (or reversibility of the scatter- 
ing paths) and the other two are enough to guarantee 
path independence for arbitrary N e m2£& With these con- 
ditions, any multiparticle scattering process can be fac- 
torized into pairwise scattering events and there is no 
ambiguity in the multiple ways of carrying out the fac- 
torization since they are all equivalent (the reader can 
find an illustration of the situation in four-particle space 
in the review article of Ref. 1^3 ). We remark that in 
the case of impurity models, any one of the indices in 
these relations can take the value '0' that stands for the 
impurity. 

These three conditions together with the constraint 
coming from the Schrodinger equation, constitute an 
over-constrained algebraic system for determining the 
electron-electron S-matrix (S = Si ,2)- Nevertheless, it 
admits a solution. 

It can be shown that the only solution is the following: 



Si, 2 — 



(fci - hi) - iV 2 V\ a (fci - fc 2 ) + iV 2 V\ 2 



[k 1 -k 2 )-iV 2 {k 1 -k 2 ) + iV 2 ■ 

This matrix serves to describe the electron band with 
linear dispersion within a basis of reduced symmetry 
SU(NM) -> SU(N) <g> SU(M), consistent with the re- 
duction in symmetry operated by the addition of the im- 
purity terms to the Hamiltonian of the host band. It is 
important to emphasize that the introduction of Si^ does 
not signify that we have modified the original Hamilto- 
nian -f/MchA by introducing electron-electron interaction. 
Instead, the choice of Si,2 corresponds to a choice of basis 
in the space of free electrons^ 



III. PERIODIC BOUNDARY CONDITIONS 
AND BETHE-ANSATZ EQUATIONS 

We proceed to impose boundary conditions. This is re- 
quired in order to be able to properly count and label the 



states. We will discuss here the case of periodic bound- 
ary conditions, imposing the following set of conditions 
on the iV e -electron wave function, 



43 ({*» 



where we are considering a finite ring of length L. As we 
are able to move electron j to the far left (xj = —L/2) 
or to the far right (xj — L/2) using the S-matrices, the 
boundary condition gives rise to the following eigenvalue 
problem: 

ZjA= ZjA 

where the eigenvalues Zj = e~ lkjL of the transfer matrix 

z i = ■ ■ • SjiSjoSjN ■ ■ ■ Sjj+i 

will allow us to find the spectrum of the Hamiltonian via 
E = kj. A is a vector in the internal space of the 
N e electron sector. Equivalently, this condition corre- 
sponds to taking particle j around the ring and asking 
that the wave function should not change; it should hold 
for all j = 1, . . . , N c . The original problem of finding 
the eigenenergies and eigenf unctions of the Hamiltonian 
is thus reduced to that of finding the vector amplitudes 
in internal space (A) that are simultaneous eigenvectors 
of all these eigenvalue problems. 

The solution of this class of eigenvalue problems was 
first tackled by Yang)S who solved the problem by 
means of a 'second Bethe- Ansatz' (see also the work of 
Baxter~). In the late 1970s the procedure was system- 
atized into what is known as the Quantum Inverse Scat- 
tering Method (QISM). 36 However, the existing technol- 
ogy is insufficient for our purpose; the structure of the 
impurity S-matrix in the multi-channel Anderson model 
requires a reformulation and extension of the standard 
formalism. We give a detailed account of those develop- 
ments in one of the appendices at the end of the article. 



Bethe-Ansatz Equations for the Multi-Channel 
Anderson Impurity Model 



The eigenvalues, zj = e 



of the transfer matrix 



are given in terms of the charge rapidities (kj), which 

in their turn, together with the spin rapidities (A„ ) 

and the quadrupolar- flavor rapidities (A„ ), completely 
specify the particular eigenstate (see first appendix) . The 

spin rapidities (A« ) and the q- flavor rapidities (A'^) 
describe, respectively, the spin and flavor dynamics as 
well as the symmetry of each state. The index x = s, (q) 

specifies that refers to a spin (q-flavor) degree of 

freedom, the index r describes the rank (related to the 
spin or flavor symmetry of the state) , and finally a labels 
the different rapidities of each type and rank. 

The charge, spin, and flavor rapidities must satisfy a 
set of equations - the Bethe-Ansatz Equations (BAE) - 



6 



that are derived in the first appendix. These equations 
encode the full information of the model: 



M i k - A s(1) - »A M? k - A q(1) + » A 
11 ■ - A s(1) + 4 A 11 fr -A 9(1) -«A 

n=l "'J il « ^ tiA m=l ""J Jt «i 4LA 

with the conditions, 



m* 

n 


A s(r) 


- A s(r) 


-i2A 


= n 

m— 1 

£T=±1 




a a(r+cr) 


-iA 




A s(r) 


- A s(r) 


-H2A 




A s(r+a) 
1 *-m 


+ iA 


n 


A 9(r) 


- A 9(r) 


-i2A 


JVf 9 , 

= n 

m— 1 
<r=±l 


A 9(r) 




-iA 




J Ml 


- A 9(r) 


+ i2A 




a g(r+ff) 


+ iA 



where for convenience we have used the definitions, 

A*, 9(0) _ jl 

K\ {M) =e 

and accordingly M°' q = N c , M q M =Ni = l, and M S N = 0. 
One sees that the effect of the impurity enters via the 
auxiliary conditions for the q-flavor rapidities. This is 
a distinguishing feature, different from what happens in 
the equations for the single-channel Anderson model, or 
in those for the Kondo model regardless of the number 
of channels. 

Each solution of the BAE corresponds to an eigenstate 
of the Hamiltonian. The eigenfunction can in principle be 
determined from the eigenvector A and the corresponding 
energy eigenvalue is given by 

£ = 5> ■ 

3 

I 



The charge-, spin- and quadrupolar-rapidities entering 
the solutions of the BAE are in general complex, and take 
the form, in the thermodynamic limit, of strings£L2&2^ 
An n-string of spin or q-flavor rapidities consists of n 
equally spaced complex numbers symmetrically arranged 
around the real axis (a; = s, q): 



A*M = A*M +i(n+l — 2a) A with a = 1, . . . , n . 



Strings are thus specified by a single real number A„ , 
in terms of which we shall rewrite the BAE. Similarly, 
there also complex values of the charge rapidities kj cor- 
responding to bound states among the bare electrons that 
build up the theory. They too form strings and will be 
incorporated in the BAE through the real part charac- 
terizing them. Consider a bound state of n particles, the 
corresponding charge rapidities form a string: 



*4 0) = fc ( "- 1} + i (n + 1 - 2a) A with a = l,...,n 

and we have spin-rapidity strings associated with them 
in the first n — 1 ranks: fc^" -1 ) = for all ranks such 
that n = r + m. 

We now incorporate the string solutions into the equa- 
tions and simplify them using the notations and relations 
given in the second appendix. After some algebra, one 
arrives at the final version of the discrete-real-BAE: 



plus, 



e 3 = 11 H e «m ( 

m— 1 i—1 




-r i] )nn^ 

m=l (3=1 


^-asoiiiici 

m=l (3=1 




n n e -» ( a - 

m=10=l 




i=i 


~ 13 ) nnn <m {kv - 

a=±l m=10=l 


.s(r+cr)\ 


n n e - ( A i 13 

m=l 0=1 


-A 9(1) " 

l ^m(3 


m—1 i—1 


- fc i m_1) ) n n < m (aj£j 

m=l (3=1 


_A 9(2)N 1 


n n e - ( A ^ r) 

m=l (3=1 


ly m0 , 


)= n nn<j 

a=±l m=l 0=1 


' A«W A«(H-<r)\ 





where we continue to use the notation A' = e. The 
set of solutions of these equations determines the set of 
eigenvalues of the problem. Now all the unknowns are 
real valued variables. 



1. Continuum Distributions Formulation 



One can study the system in the thermodynamic limit 
where its size, L, and the number of electrons, N Q , both 
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tend to infinity with the density N c /L being held fixed 
(one could later consider the scaling limit, when the den- 
sity is allowed to go to infinity while some physical scale 
is kept fixed) 42i In that case the separation of distinct 
"solutions" is of order 0(1/N C ), so that rather than con- 
sidering explicit "solutions" one can describe the system 
in terms of "densities of solutions" , typically to be de- 
noted by p(z), describing the number of solutions falling 
in a particular rapidity interval (z, z + dz)H 

To determine the densities one proceeds by taking log- 
arithms of the BAE obtaining transcendental equations 
for the rapidities characterized by integers that arise from 
the logarithmic branches. These integers label the "so- 
lutions" and are the quantum numbers of the eigenstate. 
One then constructs the counting functions^ for the dif- 
ferent rapidities [in our case v cn (z), vin, (2), and v q n (z)]. 
The counting functions range over all integers: those that 
have been selected for a state correspond to "solutions" , 
or "roots", and those integers that are omitted corre- 
spond to "holes" . We denote by p r and p h the various 
densities of roots and holes. For example pqn (z) denotes 
the density of roots of rank-r n-strings of quadrupolar 
rapidities. These density distributions are related to the 
counting functions: 

L~ l d z v cn (z) = p cn (z) = p r cn (2) + p h cn (z) 
L-^J (z) = p& (z) = p^ (z) + p h s ^ (z) 
L^d.vft (z) = p« (z) = p%) (z) + p$) (z) 

and should also obey the following relations, 

f T N cn = J kP r cn (k) 

£M&> = (A) 

. iM« = J x p#> (A) . 
These quantities are combined to define further ones, 
r N n 



E nN e 

71=1 



M ( s r) = E n M t 

71=1 



(r) 



71=1 



where N c = N e corresponds to the number of electrons 
in the system. Further, the energy (that coincides with 
the momentum for a system with linear dispersion) and 
the number of particles with different spins or q-flavors 
are given by 



T E =T^+ E nj k kp r cn (k) 

71=1 



m* r = Afi r 1} - M s (r) 



E n c . 



where r = 1, . . . , X for SU(X) (X = M, N) and we re- 
mind the reader that M q M ^ = N } = 1 is the number of 
impurities in the system. Consistently with this conven- 
tion we have the density p r q i (A) = L~ 1 5 (A — e). We 
will use these quantities to couple to crystal fields when 
we compute the thermodynamics. 

Starting from the derivatives of the logarithm of the 
BAE and using the density distributions defined just 
above plus the convolution kernels defined in the sec- 
ond appendix one can write the continuum version of the 
BAE: 

n h —Hi n * n r _ is * /i r (") 4- Ft * n r W 

Pen ~ lyZ °7im * Pern * P sm T O n m * P qm 

h(r) _ _ a * „r(r) , ts r 1 

T Dnm * Psm + a nm * P sm 

h{r) _ _ a r(r) , r r> r , 

Pqn — A nrn * P qm T O^l-Oram * PcmT 



(repeated indices are contracted). The way they are writ- 
ten, these equations determine the densities of holes as 
a function of the densities of roots and their solutions 
correspond to the different eigenvalues of the system. 

A detailed analysis of these equations will be presented 
elsewhere. In the following we will use them as a starting 
point to write a second set of equations whose solution 
allows one to compute the free energy of the system and 
gives that way access to all thermodynamic quantities. 



IV. THERMODYNAMIC BETHE-ANSATZ 

The extension of the Bethe-Ansatz formalism to ob- 
tain finite temperature information was first done in the 
case of the Bose gas by C. N. and C. P. Yang4i It was 
later adapted to the study of spin chains in the works of 
Gaudin^Sii^ Takahashi,^^ and others. The formalism 
is by now well developed and goes under the name of 
Thermodynamic Bethe-Ansatz (TBA). In the context of 
impurity models it was extensively used to find the impu- 
rity contributions to different thermodynamic quantities. 
In the following we outline the main steps and results of 
the TBA procedure as it applies to the multi-channel 
Anderson impurity model. 

We shall proceed in a standard manner and derive an 
expression for the free energy of the system. We will work 
in the grand canonical ensemble, where the free energy 
is defined &si& 



F = E-TS- fi c N c 



/ h sr TYl sr 



^ ^ hq r TYlq r 



m qr = M q 



{r - 1] Mf } + 5 ril N c 



Out of all the elements that enter this expression, the 
only one not mentioned in the previous section, and 
requiring special attention here, is the entropy. We 
define the entropy of a density distribution as follows 
(subindices are suppressed): 

±S{p r ,p h }=p r ln(l + r,)+p h ln(l + fj) 
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where we introduced the notational conventions rj = 
p h /p r and fj = rj^ 1 = p r /p h - S {p r , p h } measures the 
number of microscopic states consistent with a given 
macroscopic (thermodynamic) state given by p r and 



ph4± p or ^ ne sa k e Q £ further compactness, we also in- 
troduce the notations f = In (1 + r/) and f = In (1 + fj). 
We write down the explicit expression of a free energy 
functional for the system: 



\t {p^} = \e q + -j-hqM + E f [n (k - fx c - h ql ) -^h si -T f cn ] p r cn - £ f T f cn p h a 



(r)_ 

STL 



E / [- - h >r) n ~ T f ™] ffl E / T I^PI 



E / [~ ( v+i 



V)n-Tf« 



n r ( r ) 

I qn 



We next seek to determine the free energy, that in the 
thermodynamic limit we are allowed to evaluate as a sad- 
dle point: we vary T with respect to the distributions p r 
and p h , subject to the constraint that they must satisfy 
the BAE. We thus obtain the so-called Thermodynamic 
Bethe-Ansatz equations: 

fc„ =n (k — /J, c - h q i) /T - 5 n >ih si /T+ 

+ G n m * f"em K m * f^ m ^ B nm * 

fW =n (^ sr - W) /T + A'„ * f cr + AZ n * fW 

=n (v - v+i) /r - * fem + * fW 



(repeated indices are contracted and the kernels are 
given in the second appendix). It is possible to refor- 
mulate the TBA equations as a set of recursions linking 
the different unknown distributions^ Such a formulation 
does not involve infinite sums and has also the virtue of 
rendering the structure of the problem more transparent. 
After a few algebraic manipulations - making extensive 
use of the recursion relations for the convolution kernels 
mentioned in the second appendix - one can re-express 
the TBA equations in what we call their recursive formu- 
lation (sometimes referred to as the Gaudin-Takahashi 
form). They readt^ 



Equations for the spin-rapidities: 

fft? = lex (<5„ 4 G * f cr + G * f s ( ^ +1 + 5 n ,iG * f^ - G * f~i„ +1) - G * f£r 13 ) 
Equations for the flavor-rapidities: 

£# = lex (-8 n < N 8 r , x G * f cn + G * f^ + i + 5 n<1 G * i^-i ~ G * f^ +1) - G * 
Equations for the charge-rapidities: 

{ cn<N = lex (g * f$ + G * i$ - G * f cn+ i - S nA G * f cn -i) 
plus the 'driving' equation, 

f ( cN =\e*(N$-KGl M > M Kf^+«G^ M Ki cm ) 

\ f cN = lcx(-MA-« G ^^f cm + «<f^f;r t) ) 

I 



where the dimensionless variables A = ir(k — p)/2A and r = 7rT/4A were introduced^ and, for convenience, we 
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defined the function lex (a;) = In (1 + exp (x)). The gen- 
eral recursion kernels (Gm' M ^) and their regularized ver- 
sions { R Gm' M ^) are discussed in the second appendix. 
The expressions for f c jv and i C M can of course be derived 
from each other, but here we write down both of them 
explicitly because both are useful in the numerical stud- 
ies. 

These equations are complemented by asymptotic con- 
ditions that are also a direct consequence of the 'non- 
recursive' TBA equations: 

K , i * f (r ) - K * f (r) > rv r = hxr+1 ~ hxr 

where h xr is either /i sr =i,...,iV or h qr =i,...,M for the spin 
or quadrupolar f-distributions, respectively. 

The structure of the recursive relations between the 
different distributions can be nicely visualized in a graph. 
In Fig. we show an example of such a graph, a 'TBA- 
diagram' for a particular realization of the multi-channel 
Anderson model. 

This representation of the equations highlights their 
connectivity and is the natural extension of the graph- 
ical representation commonly used in the case of other 
integrable models (for comparison and an example of 



another related impurity model see the work on the 
multi-channel Coqblin-SchricfTcr model4&). The differ- 
ent nodes correspond to the different distributions: (i) 
those in the horizontal stripe extending towards the left 
represent the rjin , (h) those in the vertical stripe ex- 
tending downwards represent the fjj$ , and (iii) those in 
the vertical column enclosed in a box correspond to the 
fj cn . The graph was drawn for the particular case of 
SU(N = 4) ® SU(M = 5) symmetry, but its structure 
is generic. The horizontal straight lines indicate that the 
nodes are two-way connected by the equations according 
to lnr^nodo = G*Tln (1 + ?7ncighbor) and the vertical wavy 
lines indicate that the two-way connections are given by 
htffoodc = ~ G * Tin (1 + jyncighbor)- The double straight 
lines highlight the special kind of connections in the case 
of f) cN . 

There is a dual interpretation of this graph in terms 
of the reciprocals of all the distributions. Within this 
new picture the different nodes are as follows: (i) those 
in the horizontal stripe extending towards the left rep- 
resent the ffsn , (ii) those in the vertical stripe extend- 
ing downwards represent the rjqn, and (iii) those in the 
vertical column enclosed in a box correspond to the 
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r\ cn . The meaning of the straight and wavy lines also 
gets interchanged: horizontal straight lines now mean 
hirfrodc = —G*T In (1 + 7?neighbor) and vertical wavy lines 
mean ln?7 no dc = G * Tin (1 + ^neighbor)- This duality 
plays a role when one considers the two different inte- 
gral valence limits, such limits will be discussed later. 

The two boxes shown in the graph enclose those nodes 
that enter the expression for the impurity contribution to 
the free energy. Starting from the free energy functional 
and evaluating it using the distributions that obey the 
TBA equations, after a certain amount of algebra, one 
arrives at the expression F = -Fbuik 



-Pimp where 



F] 



bulk 



L 

2 



- Vn / Tln(l + fj c 



and 



Vq 



^G| M ' M) *Tln( 



1 



(*) 



(M+m,M) 



Tln(l + fj c 



from which all the different thermodynamic quantities of 
interest can be derived. Here Fbuik is the bulk contribu- 
tion to the free energy (an extensive part that recovers 
the standard result for free electrons and is all there is 
in the absence of the impurity), on the other hand -Fl mp 
is the extra contribution due to the presence of the im- 
purity. In the following we shall pursue the study of the 
latter. 

As expected for an infinite flat-band, the bulk part of 
the free energy is found to be divergent. To study it, 
some form of regularization should be introduced, e.g. a 
bandwidth cut-off. A convenient procedure that was suc- 
cessfully applied in the Bethe-Ansatz study of the single- 
channel degenerate Anderson model is the introduction 
of a 'Lorentzian cut-off scheme' J^iiSiil In that frame- 
work it is usual to fix the density and let the chemical 
potential be determined as an implicit function of the 
cut-off parameter. As a result the impurity thermody- 
namics is found to be naturally expressible in terms of 
a certain scaling invariant parameter, e* (e, A)j3^ an im- 
plicit function of cut-off, temperature, and fields4£ On 
the other hand, the impurity contribution to the free en- 
ergy is found to be regular (i.e. finite even in the limit of 
infinite cut-off) . Since we shall be interested solely in the 
impurity thermodynamics, we will introduce no cut-off 
and continue to work always in the grand canonical en- 
semble, keeping fixed the chemical potential rather than 
the density. With this convention, we can continue to 
describe the physics in terms of the original microscopic 
parameter e. 



V. IMPURITY THERMODYNAMICS 

In this section we study the impurity contributions to 
the different thermodynamic quantities of interest. We 



will first list some analytical results and then give an ex- 
tended discussion of the numerical solution of the TBA 
equations and the results obtained for several thermo- 
dynamic quantities across the different regimes of the 
model. 



A. Analytical Results 

As already pointed out in the two-channel case^ the 
multi-channel Anderson model displays a non-trivial zero 
temperature limit for the impurity contribution to the 
system entropy. This is a clear indication of the non- 
Fermi liquid character of the ground state. We can 
find the value of this entropy and identify the relevant 
scale for the cross-over into the low temperature phase in 
closed analytical form and for the general multi-channel 
case. 



1. Zero Temperature Solution 

We begin the study of the TBA equations by taking 
the zero temperature limit of the equations in order to 
identify the ground state. This is a required preliminary 
step before abording the study of the impurity contribu- 
tion to the residual entropy. 

It is convenient to introduce the distributions £ = 
Tin rj. Assuming that the derivatives cV£ are bounded 
distributions, or have at most isolated logarithmic diver- 
gencies, we derive the following limits^ 



lim Tin (1 + n) 



lim Tlnfl + r?) = f 

T^0+ 1 1 



These limits are the key step in the derivation; using 
them it is immediate to write down the zero temperature 
limit of the equations, 



NMz =C 



Mt< 



JM-t)+ r ^ 



=G- 



■(i)+ 

qn 



£,cn <N 

+ G* 



Jn)+ _ t _| 



-1 ^'M^cn-l 



?sn+l ~r On,l? sn _i 



4n = - $n<N Sr,lG * + 



G 



where the variable is the shifted rapidity z = k — (a, and 
we use the notation <5„ i = 1 — <5 n ,i. These equations 
can be solved exactly when M = N and there are no 
external applied fields. Let us start discussing the two- 
channel case (M — N — 2). Since there is only one 
rank we drop the rank superindex during this discussion. 
It is easy to see from the equations that £~> = Vn 
and £ll = Vrt > 2. It is natural to assume in the 
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zero fields case that £~ 



0, which implies £ 



0. 



Notice that after these considerations the equations for 
the spin rapidities are decoupled from the rest. We solve 
them first, taking into account the asymptotic condition 
i^n+i*^sn — K n *£ an +i — > — 2h a = (we use the standard 
definition for the magnetic field). The solution is very 
simple: put all £ sn 's to zero. We see also that in analogy 
to the spin case we can solve for the q-flavor sector simply 
by taking all £<J^'s with n ^ 2 to be zero. Only three non- 
zero distributions remain and are given directly by the 
equations 



u = c 2 +c 



2z 



U = &=G*t 



c:2 



?2 



^ql 



r2 



Thus, in the absence of external fields, the ground state 
is built out of a sea of charge-spin strings filled up to the 
chemical potential (i.e. k = fi) and a completely filled sea 
of q-flavor 2-strings. We can derive the zero temperature 
impurity valence, or level occupancy, in a closed form 
(we use the notations A = ttz/2A and J = 2 A/ (e — fi), 
notice the charge susceptibility can be easily obtained 
taking a derivative): 



4(A-tt/J) 



(A-7r/J) 2 + 7r 5 



£ 92 (A) dX 



The occupancy is integral (i.e. nP c w 0, 1) for \e — fi\ » 
A, and non-integral otherwise. It is a simple exercise 
to show that the impurity contribution to the residual 
entropy is S® mp — In \[2 - a result that is consistent with 
what was found in the integral valence limit given by 
the two-channel Kondo model. 16 We will give below a 
derivation of S° for the general multi-channel case. 

Let us now return to the case with arbitrary values of 
M and N, but always in the absence of applied exter- 
nal fields. The equations are more involved and one is 
not able to make general statements about the positiv- 
ity of the different distributions the way we did in the 
two-channel case. We proceed by making the following 
educated guess: 



1- £sn = Vn, r (this is to ensure a paramagnetic 
ground state), and it implies £,tr <N — 0- 

2- Q[n = Vn 7^ N and for n = N the distributions 
are negative, i.e. = 0. 



The equations then take the form, 
NMz = Cmm * Cn ~ °nn * 



The first equation is decoupled from the rest and one 
should solve it to determine £ c at that will later play the 
role of a driving term on the other equations. A solu- 
tion that is exact when N = M, but approximate oth- 
erwise, is given by £ c n — z [M 9 (— z) + N (z)]. When 
approximate, this solution is not accurate in the rapidity 
interval \z\ < A. This region corresponds, in the free 
energy, to the mixed valence regime that interpolates be- 
tween two different multi-channel Coqblin-Schrieffer lim- 
its (this will be discussed more in detail in Sees. IV A 31 
and I VI|) . When N ^ M, the above solution can be used 
as the initial guess for an iterative numerical scheme; 
this idea will also be useful for the numerical solution of 
the finite temperature case. The remaining two sets of 
equations can then be 'unnested' to obtain explicit ex- 
pressions for the different remaining distributions, all of 
which are fully determined by £ c at (in complete analogy 
to the two-channel case). 



2. Residual Impurity Entropy 

Regardless of its precise form, the solution for £ c jv will 
have both positive and negative non-zero parts (as is clear 
from considering the large \z\ asymptotic limits). Thus, 
the zero temperature and fields solution discussed earlier 
indicates that ^ n<N and are non-zero everywhere^ 6 - 



(r) 

In turn this means that rj cn<N and ry qI ^ are zero in the 
zero temperature limit. At the level of the TBA equa- 
tions, this has the effect of isolating the equations for the 



(r) 
Vqn<N 



distributions: 



k»7#< w = G * [Kn-1 In (l + vS+i) + 

-in(i + ^; +1 )) -in(i + ^- 1 : 



The reader should notice our earlier result, that the re- 

(r) 

lated distributions £, qn<N are all zero at zero temperature, 

merely indicates that the respective rj q n <N are finite (i.e. 
neither zero nor divergent). The equations above con- 
tain no driving terms, therefore the distributions should 
be constant functions of z. In such a case we can make 
the replacement G — > 5 (z) /2 and performing the convo- 
lutions obtain a purely algebraic set of equations, 



(V r) 



1 



(r) 



■00 



(r) 
Vqn- 



Vqn 



1 + 



that, is easy to verify, admits the solution 
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Plugging this result in the expression for the free energy 
we read off the residual impurity entropy: 



In- 



sin ZLM sin * N 

>llL N+M _ , bm N+M 



N+M 



N+M 



This formula agrees with the one that was derived in 
the integral valence limilssi using the multi-channel ver- 
sion of the Coqblin-Schrieffer model4iL£ii That model is 
the Schricffcr- Wolff limit^i of the multi-channel Ander- 
son model that we are studying. We now turn to the 
discussion of such limit. 



3. Integral Valence or Schmeffer- Wolff Limit 

In the standard approach - at the level of the Hamil- 
tonian - the Schrieffer- Wolff (SW) limit is defined as a 
truncated unitary transformation that eliminates the di- 
rect hybridization term and traces over (discards) the 
less favorable states of valence from the impurity Hilbert 
space. This limit leads to a projection onto a Hilbert sub- 
space with the impurity site permanently occupied by a 
local moment. Starting from i?MchAi one obtains this way 
the multi-channel Coqblin-Schrieffer (CS) model (the N- 
channel SU(M) model or the Af-channel SU(N) one de- 
pending on the sign of e, see below). In the following we 
will discuss how this effective limiting procedure can be 
carried out at the level of the TBA equations. 

Let us point out that the SW-transformation can be 
viewed as a first step in a more detailed Renormalization 
Group analysis (which we will not carry out in here). 
The multi-channel CS-model is therefore, when away 
from mixed valence, the naive low temperature effective 
Hamiltonian of the full multi-channel Anderson model. 
This effective Hamiltonian is relevant for the description 
of the low energy dynamics for energy scales E <C 7h • In 
other words, 



H 



MchA 



0\ — corrections in all sectors 



(with x = s,q for the magnetic or quadrupolar sector, 
respectively). In the limit of |e| — > +00, Th — * +00 
and the projection is exact (albeit with a vanishing ex- 
change constant). For any finite Th, the effects of the 
corrections, though 'subleading' to those kept in Hcs>> 
may be important at low enough temperatures if they 
combine with other operators to provide more singular 
contributions. This will be the case, for instance, when 
the degeneracy of the lower-energy impurity configura- 
tion exceeds that of the (degenerate) higher-energy one, 
as we shall stress in Sec. I VII We conclude therefore that 
the SW-transformation needs to be used with extreme 
care in circumstances like the one described above. 

We now turn to studying the integral valence limit via 
the Bethe-Ansatz. The discussion of the integral valence 
limit turns out to be simpler using an 'inverted' set of 
TBA equations (trading back the recursions on the rank 



indices for finite sums). In this alternative formulation, 
the equations read (see the appendix for the kernel defi- 
nitions), 



-f W =R r jf * 



F (») 



-&> = - s n , N t-'ae: 



+R M * 



V»+l ' ° n > lL qn-l~ 



-f — — T~ 1 AE n - 



-G 1 * - 5 n< pf8 s if c 



f(l) , f (m) n -l 7 



where we defined the effective driving terms (we remind 
the reader that r = 7rT/4A), 



aei/t 

ae:/ t 



R M *l 



cN 



R 



N 



* I C JV 



We need to discuss separately the quadrupolar and the 
magnetic scenarios. The calculations are very similar in 
both cases. We turn our attention to the former. In 
the quadrupolar SW limit {e q < e s ), 77 the SU{N) ® 
SU(M) Anderson model maps into an A^-channel CS- 
model for a local SU(M) quadrupolar impurity. The key 
to understand this limit is the study of the driving terms. 

At low temperatures - in a sense that will become pre- 
cise below - and for |A| ^> 1 (i.e. away from the intrinsic 
mixed- valence region) we can make the following approx- 
imation - motivated by our zero temperature solution for 
the ground state - that captures the leading functional 
dependence in temperature and rapidity of the distribu- 
tion associated to the maximal charge-spin bound states: 
f cA r lex(A/A/r) 0(-A) +lex(N\/r) 0(A). The ef- 
fective driving terms can then be approximated as (for 
A 3> 1 in the quadrupolar case), 



ae: ~ -r 



M 



M l 7T (M — 

sin 

A>1 47T M 



[MX]~ — > — 



r ) 2 \ 



ae: 



R^" 1 * [NX]' 



nX . 



As the temperature is lowered, both driving terms di- 
verge. But since we have that AE™ ^> AE^ due to the 
different dependence in rapidity, AE™ /r will drive to zero 
the distributions for the charge rapidities (f C ri< N ) foster 
than AEq/t will drive the f^'s. That way the fg^'s are 
effectively cut away from the other distributions (in this 

(r) 

case the fg n s) and they alone determine the impurity 
thermodynamics. The free energy is given by 



-^imp / T 



-E 



IJ-q 



2tt 



sinf 



cos $ + cosh ^ q 



where we performed the change of variables M£ = 
2ir I J — 2A and redefined the s as functions of the 
new variable £. The goal of this change of variables is to 
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remove the coupling constant (J) dependence from the 
expression of the free energy and move it into the TBA 
equations. The only explicit dependence on J will be 
in the driving terms and the same is true for the ex- 
plicit temperature dependence; this fact will allow us to 
identify the natural temperature scales of the system. 
Remark that, in terms of the new variables, the main 
contribution to the free energy comes from the value of 
the distributions around ^~0. 

a. High Temperature Scale: Consider the driving 
terms AE r s L . Changing variables we obtain, 



AE"/t — ► nA/r = 
A>1 



nir 

T 



nM 



We shall thus define the Schottky temperature scale (the 
name will find its motivation later, with the discussion 
of the specific heat): 



rs 



Ntt 
~2J 



To = 



NA 



For t < rs (and J < 1), the valence fluctuations are 
quenched and the model goes into a regime where an 
effective description in terms of a CS-model becomes ap- 
propriate. The distributions associated to the charge ra- 
pidities go to zero and isolate the ones associated to the 
quadrupolar rapidities that form a system of TBA equa- 
tions identical to the one of the CS-model. Notice that, 
as the temperature is lowered, these driving terms di- 
verge faster than the other ones that have to overcome 
a decaying exponential in the numerator (cf. with the 
caveats about the SW-transformation discussed above). 

b. Low Temperature Scale: Consider now the driv- 
ing terms AE^. Changing variables once more we obtain 

AE r /T ^l sin *( M - r ) e Z-Hr/r,) 

ql a»i N M 

where 



NM 2 _ 2JL _ 
tk = e MJ 



M 

T K = NA [ — | < —• 



We have chosen to leave a factor of 2/N outside of the 
definition of the Kondo scale in order to have a com- 
plete resemblance between the resulting TBA equations 
and those for the scaling limit of a multi-channel CS- 
model^ (that CS-model, together with a cut-off prescrip- 
tion D c ff = NM 2 A/Att 2 , is therefore the appropriate low 
energy effective theory). Since the equations match, all 
the analysis done for that model (finding the leading ex- 
ponents of the specific heat coefficient and susceptibili- 
ties, etc.) applies in this limit and we do not need to 
repeat those considerations here4& 



c. The two Scales: For the magnetic moment limit 
(e s -C s q ) we would find again two scales with the roles 
of M and N interchanged (we omit the details since the 
considerations are very similar). As a shorthand, we can 
extrapolate the two low energy scales (quadrupolar and 
magnetic) into the high energy scales of the other regime 
(for the opposite sign of e). This serves the extra pur- 
pose of providing an ad hoc interpolation between both 
regimes and across the intrinsic mixed valence region. We 
can define thus the two scales^ 



T q ee — Tq « — (M/2) 2 In l + -e-" 

TT TT \ ZTT 



2A 2A , , n2 , / M 2. 
T s ee —t s » — (A/2) 2 In 1+ — 

TT TT \ ZTT 



The role of these two temperature scales was illustrated 
in the schematic picture of Fig. [21 In terms of them we 
can further write, 



T H (e) = max {T s (e) , T q (e)} — - T s (s) 

|e|»A 



T L ( £ ) = min{T s ( £ ),r g (e)} — T K (e) 

|e|>A 



that serve as boundaries among the high-, intermediate-, 
and low-temperature regimes. At high temperatures the 
impurity is in a mixed valence state. As the tempera- 
ture is lowered, the system crosses the first temperature 
scale indicated as Th, for |e| > A this coincides with 
Ts (the chemical potential is taken to be zero to lighten 
the notation). At this point the system enters a regime 
that can be approximately described with the respective 
CS-model. At first, the system is in the 'unscreened' lo- 
cal moment regime, extending between Th and the lower 
scale Tj^ — Tk- The larger the value of |e| the wider the 
temperature window of this regime. For |e| < A the two 
energy scales merge, indicating that there is no moment 
formation. As the temperature is lowered beyond Tl, the 
moment screening commences and the physics is asymp- 
totically governed by a line of fixed points parametrized 
by the value of e/A (or some observable that depends 
on it and varies along the line, as for instance the charge 
valence of the impurity) . The different points in the line 
share the same value for the impurity entropy, and the 
same set of leading exponents of the specific heat coeffi- 
cient (i.e. 7 = Ci m p/T) or the different susceptibilities. 
However, the prefactors of the different leading terms will 
in general vary along the line and could in principle be de- 
termined by direct measurement. Thus a multi-channel 
Kondo effect takes place for any e as the temperature is 
lowered. It is amusing to note that for e = 0, in particu- 
lar, the Kondo effect takes place without moment forma- 
tion. 
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4- Asymptotic values 

In this section we will be interested in finding the val- 
ues that the different distributions (in this case the 77's) 
take when A — > ±00. This corresponds to the limit of in- 
finite |e — n\ at finite temperature. These results will be 
required in the next section, where the numerical solu- 
tion of the TBA equations is discussed. Remark that, due 
to the behavior of the two energy scales derived above, 
in the limits considered in this section one always has 
T k <g; T <g; T s and the A -> ±00 limits shall be, re- 
spectively, identified with the infinite temperature limits 
of the effective quadrupolar and magnetic multi-channel 
CS-models. 

From the results above and inspection of the equations, 
the reader can convince himself that, 



lim fc/v 

A — >+oc 



-CO 



lim f c jv = 

A — >+oo 



lim i cN = 

A — ► — 00 



lim f c jv = +00 

A — ► — 00 



We shall use this result to analyze the TBA equations 
in the asymptotic limit. The two limits, A — > ±00, are 
different and we shall consider them separately. 



tributions. We want to find the values 77 a 



a. Right Asymptotics: Let us first study the limit of 
A — * +00 that we call the right asymptotics of the dis- 

lxn = , lim Vxn 

A — > + OG 

that the different distributions take in the limit (with 
x = s,q). Except for r] c N that is unbounded, all the 
other 77's acquire finite values and for A very large can 
be taken as constants. From the limit i C N — ► +00 it 

A — >+oo 

follows that all the fj cn 's go to zero and the equations 

Mr) < Mr) 

for the 77 sn s and those for the 77 qn 's form two identical 

sets decoupled from each other and obeying the following 

recurrence relation: 



Mr) , 

Vxn+1 + 1 = I 1 



_^(r+l) 
*? xn 



Mr) 
Vxn 



1 JL ^ (r) 
1 + 1™- 



-(r-1) 
V xn 



In the case of uniform or 'Zeeman' splitting, all the a x 's 
entering the large n asymptotic condition are equal and 
an analytic solution for the asymptotic values is known. 
It is easy to show with some algebra that the solutions 



are: 



16.46,52.53 



Mr) 
V xn 



sinh ((n + r) a x ) sinh ((n + N x — r) a x ) 
sinh (ra x ) sinh ((N x ~ r)a x ) 
(n + r) (n + N x — r) 
-0 r (N x - r) 



where N s — N , N q — M and a x lost its dependence on 
the rank index. A closed analytic solution is not known 
for the case of a more general crystal field splitting (more 
general splittings might be relevant to make the connec- 
tion with the experimental systems that motivated the 
model, this point will be discussed further in the next 
section). Notice that when the values of the 77 's are re- 
quired for a particular splitting, they can always be found 
numerically. 

b. Left Asymptotics: We study now the opposite 
limit of A — * —00, in order to find the asymptotic values 

V xn — um Vxn ■ Except for fj C N that is unbounded, all 

A — > — 00 

the other 77's acquire finite values and for A very large can 
be taken as constants. From the limit f C N — ► +00 it 

A — > — 00 

<-(r) 

follows that all the ijq^'s go to zero and the equations 
for the remaining distributions form two identical sets de- 
coupled from each other and similar to those of the right 
asymptotic limit. The solution of that case can be ap- 
plied to this one if the following identifications are made 
(cf . with the discussion of Fig. |2| : 



Mr) 


Mr) 


V sn 


VsM+ 




Mn) 


V cn < N 


VsM 


Mr) 


Mn) 


V qn <N — 


VsM- 



Mr) 
VqN 

Mr) 
V qn >N 



Mr) 
Vqn-N 



These are the same identifications that would be required 
for an explicit discussion of the Schrieffer- Wolff limit in 
the local magnetic moment regime, (e s <C £ q ), in order 
to match the resulting TBA equations to those of the 
corresponding magnetic multi-channel CS-model. 

Once the right and left asymptotic limits are both 
known, one is then ready to tackle the problem of solving 
the TBA equations numerically; as we proceed to explain 
in the latter parts of this section. 



B. Numerical Solution 

In order to access the thermodynamics of the system 
at all temperatures, we shall resort to the numerical solu- 
tion of the TBA equations. This task was carried out in 
the past for other integrable impurity models 
The equations of the multi-channel Anderson impurity 
model have many similarities with those considered pre- 
viously for other types of impurities, but they present 
also important differences. We give below a brief outline 
of the numerical procedure that we developed followed 
by the results we thus obtained. 
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1. Outline of the Numerical Procedure 

We solve the TBA equations using a standard itera- 
tive scheme first introduced in the work of Rajan^4 The 
idea is to use the TBA equations in their recursive formu- 
lation, start with some educated initial guess, and iter- 
ate them until certain convergence criterion is met (this 
general scheme is sometimes known as Kepler's method; 
cf. Ref. |58j\ Our particular implementation borrows 
ideas mainly from the previous work by Costi and Zarand 
for the anisotropic Kondo model^ In that work, the clo- 
sure of the infinite set of TBA equations into a finite set 
brings in a great simplification. This is not, however, the 
case in general and we need to address the issue of trunca- 
tion of the infinite hierarchy of recursions. The standard 
procedure is to define some boundary levels (in our case 

r in 3 +i anal fgn +v f° r su hably large values n s and n q of 
the level index) and fix them. It was done in the past by 
taking those boundary distributions as constants equal to 
the average of their right and left asymptotic values that 
are known analytically. 54,58 This approximation is good 
when n s and n q are fairly large (depending on the partic- 
ular model). In the case at hand the TBA equations are 
more complicated than those for the Kondo model (in a 
way that will become clear below) and we cannot afford 
the computational cost of taking too large values for n s 
and n q . The alternative we found was to start with some 

(r) 

educated guess for the distributions up to f^„ x+1 (with 
x = s,q) that interpolates smoothly between the right 
and left asymptotic values (discussed above) respecting 
the monotonicity properties of the solution. Then we it- 
erate the equations to find new values for the different 
distributions up to the levels n — n Sl n q . At the end 
of each iteration we update the value of the boundary 
distributions [n = n s + l,n q + 1) using the fine's suit- 
ably shifted and rescaled to match the analytically known 
asymptotic values corresponding to the boundary levels. 
This procedure is repeated until stable convergence is 
achieved on the free energy function that we recalculate 
at each iteration. 

The rapidity dependence of the distributions in the 
case of the Kondo TBA equations is associated to vari- 
ations in temperature, and the two asymptotic limits 
correspond to the zero and infinite temperature limits. 34 
This is in contradistinction to the case of the TBA equa- 
tions for the multi-channel Anderson model, for which 
the rapidity dependence is related to variations of the 
coupling (~ 1/ J) ■ We will therefore require modest com- 
putational effort to compute the free energy for differ- 
ent values of the coupling and determine quantities like 
charge susceptibilities, but temperature dependence will 
require independent runs for each value of temperature 
required. Since determining temperature dependence is 
essential, our computational task becomes typically two 
to three orders of magnitude larger than for Kondo im- 
purities; depending on the range and number of temper- 
ature points desired and not counting the inherent extra 
complexity of our TBA equations. 



Once we have a finite number of equations involving 
a finite number of continuous distributions, we need to 
discretize those distributions. This is conceptually done 
in two steps, the first one being the introduction of cut- 
offs on the rapidity axis. We have to choose large right 
and left cut-offs that enter once the distributions are ap- 
proximately constant functions of the rapidity reaching 
their respective right and left asymptotic values. Second, 
we need to discretize the axis interval between the cut- 
offs. This is done defining three regions. First a small 
region centered around zero that we discretize using a 
fine mesh. The size of this region is chosen depending on 
the parameters of the problem so that it encloses all the 
intervals where different distributions show rapid varia- 
tions. This typically happens for region boundary values 
of the rapidities such that the magnitude of their asso- 
ciated coupling corresponds to a Kondo temperature of 
the order of the temperature set for the system. Second, 
the two regions to the right and left of the central one 
are discrctized logarithmically until reaching the cut-offs 
defined above. In these shoulder regions all the distribu- 
tions vary slowly as they attain their asymptotic values. 

The different convolution kernels required in the cal- 
culations should be evaluated, this is done a single time 
at the beginning of the calculations when they are stored 
as matrices. The kernels are relatively rapidly varying 
functions of the rapidities as compared with the TBA 
distributions, so we discretize them using a denser mesh. 
We use a mesh that is locally n times (typically n ~ 10) 
finer than the one used for the TBA distributions and 
these are interpolated to the points in the finer mesh us- 
ing cubic splines. The extra points in the finer mesh are 
chosen to be the n-th order quadrature points of each 
subinterval of the coarser one. The convolution integrals 
can thus be split into the different subintervals and car- 
ried out with high precision in each one of them using 
Gaussian quadrature (corrections for the introduction of 
the right and left cut-offs are also implemented). Since 
accuracy is very important for certain applications^ the 
evaluation of the different kernels requires a high preci- 
sion implementation of the sums involving four digamma 
functions of complex arguments (see appendix) . To carry 
it out we wrote an algorithm based on 'a precision ap- 
proximation of the gamma function' due to Lanczos^ 

With everything set to start the iterative evaluation 
of the TBA equations, one last point is worth mention- 
ing. The evaluation of the right-hand side of these equa- 
tions requires the implementation of the ubiquitous func- 
tion lex (a;). This implementation is, however, unsur- 
prisingly subtle. We first notice that, since lex (x) — 
[x] +lex (— |a;|), it suffices to implement the function for 
negative arguments. This we do by first evaluating e - ' 1 ' 
and then using a careful implementation of In (1 + x) 
with the property of cancellation of rounding errors. 60 
This algorithm ensures monotonicity, machine precision 
accuracy, and 'graceful underflow'. 

We show below some of the results that were obtained 
with different implementations of the described numeri- 
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FIG. 4: Impurity entropy as a function of temperature. The different curves correspond to: e/A = (dashed lines); e/A = ±2, 
±4, ±6 (light lines); and e/A = ±8 (dark solid lines). The different panels give the results for different values of N and M 
as indicated to the right and above, respectively. Curves for different signs of e are degenerate in the diagonal panels. In the 
off-diagonal ones they can be identified from the value of the entropy in the intermediate plateaux (see text). 



cal procedure. 



2. Numerical Results 

As stated above, one of the important issues that 
makes the numerical analysis of the TBA equations re- 
quired is the study of the temperature dependence in its 
full range. We give below illustrative results of this de- 
pendence where the presence of the two cross-over scales 
discussed already on analytic grounds is clearly observed. 
We start by showing results for the model in the absence 
of crystalline or external applied fields. 

In Fig. Q] we show several plots of the entropy as a func- 
tion of temperature. The different panels correspond to 
different symmetries of the model; N = 2, 3, 4 row-wise 
and M = 2,3,4 column- wise as indicated. The different 
curves in each panel correspond to different values of the 
energy splitting between the two impurity configurations 
(e = e s — e q ) and, without loss of generality, the total 
chemical potential is taken to be zeroS In all the pan- 



els, the dark solid lines correspond to e/A = ±8 and the 
dark dashed line corresponds to the extreme mixed va- 
lence case of e = (or e = fi). The remaining thin lines 
correspond to intermediate values of the energy splitting 
(e/A = ±2, ±4, ±6) and are given to illustrate how the 
system interpolates among the different limits of large, 
positive and negative, and vanishing e. On the one hand, 
in the high-temperature limit, the entropy is in each case 
given by S^ p (N, M) = fee In (N + M); as expected since 
N + M is the total size of the impurity Hilbert space. On 
the other hand, in the low-temperature limit, the value 
of the entropy tends to the S-* mp (N, M) values we found 
analytically (see Sec. IV A~2l . We remark that both these 
limiting values are independent of e. 

For intermediate temperatures, the figure shows how 
the impurity entropy is quenched from its high- to its low- 
temperature values as the temperature decreases. This 
happens as a two stage process for large values of e (solid 
lines) and as a single stage process in the mixed valence 
case (dashed lines), all in accordance with the theoreti- 
cal discussion given above (see Fig. and the discussion 
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FIG. 5: Impurity specific heat as a function of temperature. The different curves correspond to: e/A = (dashed fines); 
e/A = ±2, ±4, ±6 (light lines); and e/A = ±8 (dark solid lines). The different panefs give the resuits for different values of N 
and M as indicated to the right and above, respectively. Curves for different signs of e can be identified from the iocation of 
the Kondo anomaly or by comparison with the entropy plots. 



in Sec. I V A 3 cfl . Following the solid lines as the tem- 
perature lowers, the first quenching stage corresponds to 
the cross-over scale Tg and the impurity entropy attains 
the values ks In (N) or fes In (M) for e s <C e q or e s S> e q 
respectively. (As expected, in the case of N = M, the 
curves for different positive and negative values of s are 
degenerate; on the other hand, for N ^ M, the curves 
for the same absolute values but opposite signs of e are 
in precise correspondence upon exchange of the values 
of N and M) iSl These intermediate- regime plateaux cor- 
respond to the formation of a free local moment when 
|e| S> A. As the temperature lowers further, follow- 
ing always the solid lines, the systems reach the Kondo 
cross-over scale, Ik, below which the entropy tends to the 
universal values Sf mp (N, M) characteristic of the differ- 
ent infrared non-Fermi liquid fixed points. This second 
quenching stage takes place as electrons in the different 
channels compete to screen the local moment of the inter- 
mediate regime. It is this dynamical frustrated screening 
process that is responsible for the non-trivial nature of 
the fixed point. In the mixed valence case, exemplified 
by the dashed line, the two cross-over energy scales are 



of comparable magnitude and as a result the quenching 
of the entropy as the temperature is lowered happens on 
a single stage. Remarkably, the same limiting value of 
the impurity entropy is found also in this case; all con- 
sistent with the picture given above, that the infrared 
physics is governed by a line of fixed points (Sec. IV A 3 cfl . 
In a previous work, we explored this explicitly (in the 
N = M = 2 case) using Boundary Conformal Field The- 
ory and showed that the different points in the line are 
connected by an exactly marginal operator m& 

The results for the impurity contribution to the spe- 
cific heat - obtained upon differentiation of the impurity 
entropy - are shown in Fig. [S] The different panels, and 
the different lines in each of them, follow the same con- 
ventions as in the entropy plots. In accordance with the 
generic two steps shape of the entropy (degenerate in 
the mixed valence case), the specific heat shows a gen- 
eral two humps structure (again degenerate for mixed 
valence). The lower temperature one is sometimes re- 
ferred to in the experimental literature as the 'Kondo 
anomaly', and its location is correspondingly given by 
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FIG. 6: Impurity charge valence as a function of temperature. The different curves correspond to: e/A = (dashed lines); 
e/A = ±2, ±4, ±6 (light lines); and e/A = ±8 (dark solid lines). The different panels give the results for different values of N 
and M as indicated to the right and above, respectively. Curves for e positive or negative are, respectively, below or above the 
e — line. 



Tk- On the other hand, the higher temperature one is 
referred to as the 'Schottky anomaly' and its position is 
given by Tg. There is of course no new information in 
this figure as compared with the previous one, but it is 
the specific heat rather than the entropy or the free en- 
ergy the quantity that is most often 'directly' accessible 
in the experiments. 61 

It is illustrative to look as well at the behavior as a 
function of temperature of the impurity charge valence, 
n c ,im P = J2a {ftU)- This is provided in Fig. El following 
identical conventions as in the plots for the entropy and 
the specific heat. As expected, in the high-temperature 
limit, the impurity valence approaches the values, 



n~ imp (JV,M) 



N 



N + M 



corresponding to the impurity Hilbert space fraction with 
'magnetic character'. This is an expression of the fact 
that, in that limit, all the impurity states are equiprob- 
ablc. Subsequently, as the temperature crosses the 



value Tg, the impurity charge valence changes and ap- 
proaches rapidly what will be its zero temperature value, 
n c,im P W M). This change is at the origin of the Schot- 
tky anomaly in the specific heat. As it should be, the 
impurity approaches integer valence when |e| 3> A. The 
charge valence goes to zero in the quadrupolar limit of 
large and positive e, and goes to one in the opposite, mag- 
netic limit, of e large but negative. The values of n® imp 
for intermediate energy splittings are difficult to calcu- 
late and would constitute a good non-trivial test for ap- 
proximate theories like those based on (l/iV)-expansions. 
Careful comparisons of NCA and NRG were carried out 
this way in the single-channel case^ but for the multi- 
channel case the NRG calculations rapidly become very 
demanding for present day computational resources and 
this kind of comparisons were not done. Also, since the 
full cross-over takes place at the highest energy scale, the 
convergence of n° ; should be relatively fast in NRG 
calculations. Thus in the future it might constitute a 
useful observable to monitor the progress of such compu- 
tations by comparing with the exact solution. 
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As a check, we verified numerically that 

< imp (N,M)\ s +nl imp (M,N)\_ e = l 

so that, in particular, one finds that n° c imp (N, N) | = 
1 /2. The analytic form of n° imp (iV, M) can be computed 
using the zero temperature results of Sec. IV A 11 (with the 
caveats given there for the case of N ^ M). For instance, 
in the two-channel case, we have 

n° (2 2)--/ +0 °/° Z dz x 

c ' imp ' I" J-oo i-oo cosh (a; - 2;)] 

— e + x 

x 2 "- T 

(^-£ + 2;) 2 + 4 A 2 ] 

{« 1 for £ -C -A 

= i for e = 
w for e > A 

(where we reintroduced the chemical potential). The in- 
ner integral is related to £ g 2 and has an exact expression 
in terms of dilogarithms. It is interesting to observe how, 
for certain values of the energy splittings, the charge va- 
lence first increases and then decreases as a function of 
the lowering temperature (or vice versa) as it goes from 
n* imp to n° imp . This behavior is already present in the 
case of single-channel degenerate Anderson impurities. 62 

Before proceeding, let us add some remarks on the 
technical side. We calculated both the entropy and the 
specific heat by numerically differentiating the free en- 
ergy computed using the TBA distributions obtained 
with the algorithm described in the section above. Nu- 
merical differentiation is a delicate procedure very sus- 
ceptible to truncation and roundoff errors and in gen- 
eral sensitive to any noise in the original data. The spe- 
cific heat is, since a second derivative is involved, rather 
sensitive to this type of errors (the discretization of the 
temperature plays an interrelated role as well). In the 
event of better accuracy being required, alternative ways 
of computing derived thermodynamic quantities are pos- 
sible. One common procedure is to set up secondary sets 
of integral equations for the different derivatives of the 
TBA distributions (see, for instance, Ref. 56]). Once the 
original distributions are found, these equations can be 
solved to determine their derivatives. Using them one can 
calculate directly (or with a smaller number of numeri- 
cal differentiations) the sought derived quantities. The 
procedure to solve these auxiliary sets of equations will 
be in general similar to that used for the original TBA 
equations; making the total computational effort increase 
accordingly. 

We turn now to discuss the effects of external fields on 
the physics of the impurity. For the sake of simplicity 



we restrict ourselves to the two-channel case (N = M = 
2). Since magnetic and quadrupolar fields are relevant 
perturbations, the presence of any of them has important 
effects on the entropy. In fact these perturbations drive 
the system to a totally different - this time Fermi liquid - 
line of fixed points characterized by a zero value of the 
residual impurity entropy (cf. with the situation in the 
single-channel case). This is illustrated in Fig. [7| where 
the entropy as a function of temperature is shown for 
different values of energy splitting between doublets (e) 
and quadrupolar field (h g ). 

In the upper left panel of the figure we reproduce 
again the results for zero field. The curves with posi- 
tive and negative s coincide in this case (remember we 
are setting fi = in this discussion). As we go from 
high to low temperatures, the system crosses over from 
a state with entropy Si mp = fes In 4 to one with entropy 
Simp = &b hi V2. This second value is one of the hall- 
marks of the non-trivial non- Fermi liquid fixed point that 
governs the low energy physics of the models In the 
next panel to the right we show the effect of turning on 
an external field. For h q /A = 10 -6 one observes that the 
positive-e curves stay in an Si mp = &b hi v2 plateau for 
a very short temperature interval (that disappears alto- 
gether when the impurity energy splitting is sufficiently 
large) and a third quenching step of the entropy takes 
it to a zero value final state, indicative of a Fermi liquid 
fixed point. As the reader can observe across the different 
panels, the quenching to zero entropy of the large posi- 
tive £ curves takes place when the temperature is lowered 
until it becomes of the same order of magnitude as the 
applied field (T ~ h q ). 

On the other hand, the fate of the negatives curves be- 
comes more evident as we increase the field further (see 
the middle two panels). In this case too, the non- Fermi 
liquid fixed point is unstable and the entropy goes to zero, 
but this takes place at much lower temperatures than for 
positive £. What happens is that at low temperatures 
the splitting of the higher multiplet (that plays the role 
of the orbital channel in the usual multi-channel Kondo 
scenario) renders the model an effective spin-half single- 
channel exchange model in which the impurity is 'ex- 
actly' screened at low temperatures and the fixed point 
is a Fermi liquid^ Curves with intermediate values of 
£ interpolate continuously between these two behaviors 
very much in the way they did in the other cases that we 
discussed above. 

As the field is increased and approaches the scale of the 
first stage of entropy quenching (i.e. T$) and beyond, the 
different-£ curves start to collapse and the temperature 
of the cross-over transition moves up following the field 
(compare the last two panels in Fig. 01 . The entropy is 
quenched in a single stage from S- lmp = In 4 to zero at 
the same time as the impurity charge valence goes from 
"■cimp = 1/2 to zero. 
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FIG. 7: Entropy as a function of temperature for various applied fields (we have chosen to consider quadrupolar flavor fields). 
The top-left panel corresponds to the zero field case. To the left of it and following downwards the field increases hundred-fold 
each time as h q = 10 _6 A, f0 2 A. The values of e are the same as in Fig. 0] (the chemical potential was reintroduced in the 
notation for this figure), and the top-left panels of both figures are in correspondence. 



VI. SUMMARY AND OUTLOOK 

We have carried out a detailed analysis of the multi- 
channel Anderson impurity model. We have demon- 
strated its integrability, discussed the details of the 
Bethe-Ansatz solution, and shown a number of illustra- 
tive results for the thermodynamics of the model. 



The Bethe-Ansatz allows the study of the model on all 
energy scales, we proceed now to summarize our findings 
and to put them into context. We begin by discussing 
the low-temperature regime and the characterization of 
the low energy fixed point theory that corresponds to the 
microscopic model. We have identified the existence of a 
line of boundary critical fixed points that governs the low 
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energy physics of the impurity (illustrated in Fig. [8J and 
shown how it is related to the fixed point theory of the 
multi-channel Coqblin-Schrieffer model. In fact, every 
point in the line describes non-Fermi liquid physics. At 
the extreme ends (|e| — > +00), the fixed point theory 
corresponds to an SUm {N) or SUn (M) Hamiltonian 
(in the sense of BCFT). For any finite value of e, and in 
particular for all the values that correspond to the mixed 
valence region, the physics is given by a line of SUm (N)(B 
SUn (M) theories. Along this line of fixed points the 
behavior, for instance, of the specific heat would be given 
by (cf. Ref. H), 



SUm(N) 
• 



H*(e/A) 



-co 



E/A 



SUn(M) 
• 



CO 



N=£M 



X 2 c (e) T + A 2 (e) T«+m + A 2 (e) T~ 




(N = M corresponding to the marginal case when both 
spin and quadrupolar sectors contribute T In T leading 
temperature dependences; cf. Ref. |30|). The smallest 
exponent dominates as T — > 0, except possibly at the 
limit when e — > ±00 and A 2 (e) or A 2 (e) vanish, re- 
spectively. Thus these two limits do not commute when 
/ < n for the SU / (n) end of the fixed points line (with 
n, f = N,M > 1). Notice therefore, that as e is var- 
ied along the line at small but finite temperature, the 
observed critical behavior will vary accordingly. 

The same finite-T cross-over in non-Fermi liquid char- 
acter (magnetic vs. quadrupolar) would manifest itself 
in a comparative study of both susceptibilities, whose 
expected leading low-temperature behaviors, after sub- 
traction of asymptotic temperature-independent contri- 
butions, are: 

r> , . N-M 

- x s o ~ A 2 e 

o , . M-N 

■ Yn ~ A 2 (e) T"+i« 

(here N = M is again the marginal case, for which 
InT dependencies are expected). Set for instance e in 
the quadrupolar regime, then the corresponding A 2 , (e) is 
large while A 2 , (e) can be arbitrarily small. However, one 
will still find that, for all finite values of e, the spin sus- 
ceptibility will eventually dominate over the quadrupolar 
one if it happens to carry the singular exponent (i.e. for 
N<M). 

Note that the level with the lower degeneracy always 
dictates the physics at sufficiently low temperatures even 
if it is a very high-energy level. This may be surprising 
from the point of view of the Schrieffer- Wolff limit (see 
Sec. IV A~3t . Naively, one would expect that the effects 
of the energetically unfavorable level could be simply in- 
tegrated out. This is not always the case. When the 
degeneracy of this level is the smaller one, it will end up 
dominating at sufficiently low T (as long as |e| remains 
finite) i£i In other words, discarding the energetically un- 
favorable yet less degenerate configuration also eliminates 
the frustration that is induced in the bulk electrons when 
they try to screen the virtual moment of such an excited 
state; but it is this frustration that would have been re- 
sponsible for the appearance of singular exponents in the 
impurity thermodynamics. 



FIG. 8: Schematic representation of the line of fixed points. 



This picture might help shed light upon some of the un- 
explained and sometimes contradictory results observed 
in the experiments^^ In particular, interesting possi- 
bilities are open up for better understanding of the 
intermediate-to-low temperature phase of certain heavy 
fermion compounds, inespecially those believed to be 
near mixed valence (the list is rather large, for an ex- 
ample see below). 

Thus far we discussed the low-temperature regime (FP 
- dominated by the line of fixed points). We have also 
studied in detail the other regimes: the high-temperature 
valence fluctuation regime (FV) and the intermediate lo- 
cal moment regimes (LMM and LQM) — the initialisms 
refer to the names used in Fig. [21 We identified the two 
energy scales (T s and T q ) associated with the spin and 
quadrupolar degrees of freedom, that cross each other in 
the intermediate valence region (|e| < A) and interchange 
roles as the high-temperature (Tk) and low-temperature 
(Tl) scales that indicate the transition zones among dif- 
ferent regimes. In the intermediate valence region the 
two scales 'coincide', indicating the direct transition be- 
tween the high- and the low-temperature behaviors. The 
system never develops a local moment, but goes directly 
into the low-temperature, multi-channel Kondo like, non- 
Fermi liquid phase governed by the corresponding fixed 
points. In all the cases, i.e. for all values of e, the fixed 
points arc unstable to the application of an external field 
acting either on the magnetic or the quadrupolar degrees 
of freedom. 

In subsequent work we intend to use the results 
presented here to analyze the different multi-channel 
scenarios for certain heavy-fermion compounds like 
Ui-zThjBeisiS&Si It was shown that the two-channel 
Anderson model is not sufficient to account for the rel- 
evant number of impurity degrees of freedom required 
to explain the available specific heat measurements. 29 
More complicated impurity models, like for instance an 
SU(2)®SU (5) model in the presence of crystal field split- 
tings, hold a considerable promise in that respect j 29 i 65 
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APPENDIX A: QUANTUM INVERSE 
SCATTERING METHOD 



where tr^ denotes taking the trace in the auxiliary space 
A, and it follows that 



In the following we will skip the details (and assume 
them known) of the standard case of models with SU (2) 
internal symmetry, and refer the reader to the presenta- 
tion given in the ICTP lectures by one of the authors^ 
We will adopt a notation very similar, though not iden- 
tical, to the one in those lectures^ We will discuss di- 
rectly the case of SU(3) symmetry conveniently general- 
ized to the situation when both 'particles' and 'antiparti- 
cles' are present. This generalization is required to solve 
the eigenvalue problem of the multi-channel Anderson 
model (the standard Nested Bethe-Ansatz (NBA) as dis- 
cussed in the literature, see for instance Ref. |6q|. does 
not suit this case), and was done before only for SU(3) 
spin chainsiSLS Finally, we will close the appendix with 
the generalization to the case of SU(N) symmetry. 



1. Models with 5(7(3) Symmetry, Impurities, and 
Periodic Boundary Conditions 

Let us consider an integrable model with SU(3) sym- 
metric scattering matrices given by, 



S 



(otj - a„) I jn + icPj 

(ay - a n ) + ic 



Sj n (otj 



a r , 



SjQ — IjO 



-i5{aj— <*o) _ 1 



T ( a j) = Z j ■ 

As the amplitudes A should be simultaneous eigenvec- 
tors of all the eigenvalue problems Zj , it is necessary and 
sufficient that [Zj,Zi] = 0. We can go further and re- 
quire [T (a) , T (/?)] = for all values of a and (3. This is 
guaranteed if there exists a matrix R such that, 

R AB (a - (3) E A (a) E B (/3) = E B (/3) E A (a) R AB (a - (3) 

we will refer to this identity as the fundamental commu- 
tation relation (FCR). 



b. Fundamental Commutation Relations 

By repeated application of the Yang-Baxter relations, 
it can be shown that the matrix 



Rab (a) = S AB (a) 



aI AB + icP AB 

a + ic 



QjO = SjQ (otj - Oq) 



satisfies the FCR. If we write down the expression for the 
monodromy matrix explicitly in the auxiliary space, we 
have (using a notation that is convenient for the Nested 
Bethe-Ansatz): 

A (a) B 2 {a) B 3 (a) 
H(a)= | C 2 (a) D 22 (a) D 23 ( 

C 3 (a) D 32 (a) D 33 (a) 



We choose to consider a case where SjQ has a different 
structure than Sj n ^Q. In particular we consider the kind 
of impurity S-matrix that arises in the flavor sector of the 
three-channel Anderson model. The parameters otj , ao 
are arbitrary at this point, and will be chosen later to 
be those that specify the multi-channel Anderson model. 
The eigenvalue problem given in terms of the above S- 
matrices is not tractable with the standard NBA formal- 
ism; we develop below the required extensions. 



T (a) = A (a) + D 22 (a) + D Ai (a) . 

Fully expanding the FCR in the two auxiliary spaces, it 
can be seen, after some algebra, that the sub-matrix D 
verifies the corresponding FCR of the SU(2) case, 

Rf B (a - /3) D A (a) D B (/?) = 

= D B (/3) D A (a) R% (a - (3) 



a. Monodromy Matrix 

As a first step we define a matrix that captures the 
monodromy conditions of our eigenvalue problem and use 
its transfer matrix to rewrite the problem: 

(a) = S AN (a - a N ) ■ ■ ■ Sai (a - a{) S A0 {a - a ) 



with 



Rab («) = S A % (a) 



(2) 



^ab T LLtr ab 



This means that Yang-Baxter is obeyed by the sub- 
matrix. Defining u a = a ~' c and v a = — , one can 
write the following commutation relations between dif- 
ferent components of the monodromy matrix: 83 



(this is equivalent to adding an auxiliary extra particle 
A = N + 1). The purpose of the auxiliary space will be 
to allow us to conveniently organize the products of the 
S-matrices (see below). We define the transfer matrix as: 



A a B% 



B%A a 



Va-[3 
(2) 



BtAr- 



D a B% = u fj - a B%D a S A ' s + v^ a B s a DpP 



>(2) 

As 



(2) 



T (a) = tTA Sa (a) 



They will be used extensively in what follows. 
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c. Reference State or Pseudo-vacuum 

We use the following local basis for the individual 
Hilbert space of each particle (where the 'anti-particle' 
one is for the impurity Hilbert Space): 



Using the FCR we find (the nomenclature of wanted and 
unwanted terms is standard, and the same as in Ref. |40jp. 



ft 



wanted 



/ M 

n 

\n=l 



a-0„ 



ft 



\d),\d) 



and we define the following highest weight reference state 
which we will use to construct the Fock Space for the 
problem: 




Using the definitions a n = ° and b n 



we can 



x+ic a a+ic 

write down explicitly the scattering matrices in auxiliary 
space and see how they act on the reference state. Using 
that information we find for the monodromy matrix: 



E(a)\co) 



\oj) * 



with 



A" (a) \w) * 
A'" (a) 



N 



A"(a) = Y[a(a- aj ) 

3=1 

A'" (a) = 1±1- 



A" (a) 



Notice, for later use, that A" (atj) = 0. We define 
A (a) — A" (a) + A"' (a). For the transfer matrix we find 
(verifying that the reference state is indeed an eigenstate 
of it), 

T a \lu) = (A a + Df + Df) \u) = (1 + A„) |w) • 



d. Descendant States 

We will construct descendant eigenstates from the 
reference eigenstate that constitutes the highest weight 
state of the largest possible representation. Let us con- 
sider a state obtained by acting with the linear combi- 
nation of M flavor-lowering operators (where X Sl ,_ SM is 
an arbitrary tensor) : 



p) = J2 BS1 (ft) • • 

{«<} 



B SM {[3 M ) X s 



(ft 



wanted 




xB* 1 .. .B% M D a S { f> . . . X S1 

pi pm " Asm Asi 



M 



ft 



wanted 



\n=l 



xB£...B™T£> X Si ... Sm \uj) 



5 /3a/ J " 

where in the last line we made the following definitions: 



S (2) (a) = D a S 



(2) 



.s 



(2) 
Asi 



T^(a)^tr A E^(a) . 
Let us define the reduced monodromy matrix, 



S( 2 > (a) = sf SM (a - ft M ) . . . SZ (a - ft,) 



(2) 



A (a) B(a) 
C{a) b\a) 



written in the reduced auxiliary space, and whose ele- 
ments act only on the space of indexes {si} (i.e. that 
spanned by all the possible X Sl ... SM ). Notice that the el- 
ements of S^ 2 ' and those of D commute with each other. 
We write down their 'combined' product explicitly: 



7(2) _ r> 5(2) 



A {2) R (2) 

r m n (2) 



Since both D and ^ 2 - ) satisfy the Yang-Baxter relations, 
we have a new set of FCR that are obeyed: 

B%l{ a -[3)Ef {a)~f {13) = 

= ~%\ft)~ { l ] {a)B}l ] B {a-ft) . 

The eigenvalue we are seeking involves more than a 
particular realization of the auxiliary tensor X S1 _.. SM (cf. 
with the NBA formalism). We can define a highest 
weight reference state for the space of {si} in the stan- 
dard way (notice that S^ 2 ) is the usual monodromy ma- 
trix that appears in the SU{2) case). We therefore con- 
sider the 'combined' reference state 



M 



with 



IT), 
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given by the direct product of a reference state in the 
space of indexes and \oj) the previously defined reference 
state. In particular we have: 



A a \u) 



with A Q = Il^i a i a ~ Pn)- We can act on |w} 2 with 
the 'combined' monodromy matrix: 



si 2) |«> s 



A Q |w) s 



A a A a \lu} 2 



Choosing the reference state \Q) is equivalent to consid- 
ering only B 2, s and no B 3, s when building eigenstates. 
This is the trivial case when the electrons do not hy- 
bridize with the impurity (we want instead to include 
non-trivial eigenstates involving some B 3, s as well). In 
the more general case that concerns us we write instead 
(inspired by the form of the unwanted terms): 



= B\ 



(2) 

Let us denote by t a the eigenvalue of T (a) and by t a 
the eigenvalue of (a) = Aa ' + Da^ . We need to solve 
the auxiliary eigenvalue subproblem: 

?1 2) It) = 4 2) It) 

(notice that both T a 2 ^ and \j) are, although it is not 
explicitly indicated, functions of (i). Using the results 
from the SU (2) case we write down the wanted terms: 



4 2) 17) 



= A" 



wanted 




wanted 



- A A'" 



M 

n 



•ip- 



It) 



and also the generic form of the unwanted ones: 

4 2) It) 



unwanted 



M' 



= A ", w «-7, ( II U 'U-lv | l«Tl • • -T? • • ■ TM') 



^i 2) It) 



unwanted 



M' 



'7p-7 ? 



l"Tl • ■ ■ TM ' ) ■ 



The cancellation of the unwanted terms gives a first set 
of auxiliary conditions: 



m' . . c m a 

n 7<? ~ 7 P - ^ = Tg - QQ +»2 TT Tg - 



7g-7p + ic Tg - + i|c ^ T« - + ic 



whereas the wanted terms give the auxiliary eigenvalue: 









[ n u «-7p 


+ A-A Q 











(notice that, since A a _ = 0, the auxiliary eigenvalues 

(2) 

vanish, i.e. t a - = 0). 

Having solved the auxiliary nested problem we can go 
back to the original eigenvalue problem that we are trying 
to solve. The combined wanted terms are: 



= Bi 1 . . . B 



wanted 



0M 



' M 
IT Ua ~P" 



\n=l 



And we find that the eigenvalues of Zj = T a read as 



M 



M 



-a I'a 



\n=l 



17) = t c 



M 



Zj — tdj — 



n^-^ = n 



n=l 



n=l 



Otj - fi n - ic 

Otj - (in 



We are only left with the task of taking care of the unwanted terms. Generic ones read: 



M 



unwanted 



A a {b%B£ ...B£...B%) . . . S^ lSn ) 17) 
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ir 



. Bl M 

Pm 



(^1--.^ 1S „)|7) 



and 



unwanted 



= D a (BfcB* 



^•••^)(^1...5£)_ 15 J|7) 




r ") few 



...5l 2 ) 



so that 



tr A D a 



unwanted 




,5-)(^l...^) lS „)|7) . 



Since = 0, we have the following expression 

(2) 

for the eigenvalue of the auxiliary problem: = 

^/5 n (ilpii M /3 n -7 P ) • Combining the two contributions 
to the same unwanted term and asking that it should 
vanish, we find a second set of auxiliary conditions: 



M 

n 



Pn Pm ^ c 

/3„ - /3m + ic 



n 

3=1 



M 

n 



e. Bethe-Ansatz Equations 



Let us collect the different equations, rearrange them 
and highlight the final result. Recalling that Zj — e~ lkjL 

and performing the standard shift (3 n = — if, we 
write the eigenvalue equation: 

Mi 



e ik 3 L = J| «J_ 

71=1 "I 



A«-ii 



J "1 "2 

and after shifting 7„ = A< 2 > — ic, we write down also the 
auxiliary conditions: 



Ml A (1) A (1) ir 

A (i) _ A (i) v~ 



N A (1) 

3=1 An 



a, + if 

M 2 . (1) _ .(2) 



n 



a(1) _ A (2) , •< 
— 1 Jin il m i t r 



and 



M 2 (2) 
11 a(2) 



A 



(2) 



zc 



A^ 2) - oo - 



A^ } + ic 



i (2) 



ct - 

Mi 

n 



A 



2 

(2) 



A 



(1) 



AJ ; A (2) - A liJ 4- i 

n=l ^Ml T * 

where we have taken Mi = M and M2 = M'. 



(i) 



2. Generalization to Models with SU(N) Symmetry 

It is straightforward to generalize these results and 
write down the Bethe-Ansatz equations for the more 
general case of SU(N) internal symmetry. We use the 
notation Mq = N c for the number of electrons and 
Mn = N[ (= 1) for the number of impurities. We de- 



fine A 



(0) 



OL n (= k n ) for the charge rapidities and 
(= e) (i.e. the impurity rapidities). Then 
we write the eigenvalue equations: 



A (N) - e 



ik n L 



= u 



Mi .(0) 



A m } +ii 



L A (0) - A (1) 

m=l ilm 

and the nested auxiliary conditions: 



M r . (r ) 



n 



A^-A 



(r) 



M r 



A 



(r) 



A r ^ ) + ic 



- n n 



(r) 



,(r) 



mjtn iln ilm I" tc ct=±1 m=l A n nm 

where the ranfc varies in the range r = 1 , . . . , N — 1. 



a. Generalization to Models with SU(N) <g> SU{M) 
Symmetry 

We start with the same auxiliary eigenvalue problem, 
but with the scattering matrices: 



(a j - a n ) I s jn - ic~P s jn (a , - a n ) P. + icP? 



(ctj — a n ) — ic 

e -i5(otj-a ) _ 1 



(ctj — a n ) +ic 



S j0 -I q j0 + ;/ 



We choose to consider a case where Sjo acts non-trivially 
in the 'channel' degrees of freedom only (i.e. in 'q-flavor' 
space). This problem amounts to taking the one dis- 
cussed above and adding an extra 'isospin' to it. The 
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monodromy matrix can be written as a direct product: 
S (a) — S s (a) (a), and the transfer matrix becomes 

T a = T°T< = (A s a + D%) {Al + Dl) . 

The different steps go through as before and we get the 
Bethe-Ansatz Equations (BAE): 



kernels: 



AkjL 



Mi 

n 

n=l a j 



otj - A 



A* (1) + i% 



n 



9(1) 



2 m=l "J 4 2 



with the conditions, 

M" 

n 
n 





- A s(r) 


— ZC 


A s(r) 


- K s[r) 


+ ic 






— ic 


A*) 


- A q( - r) 


+ ic 



M r+c » s(r) a s(r+<x) _ . < 

11 X s(r) _ ,s{r+a) .< 

CT =±1 m=l JV ™ iim "i 4 s 

f CT / 

r _ 



■A/'' 



n n 

er=±l m- 



A„ W - A 



q(r+a) 



— l-c 



- A "\ r ) _ A 9 
1 Jin Ji-r 



where for convenience we have used the definitions 
A 1 — ag 

and accordingly M S ' 9 = N c , M q M = = I, and = 0. 
One sees that the effect of the impurity enters via the 
auxiliary conditions for the q- flavor-rapidities. 

The solution thus far was general. To specify it to our 
model we take ctj — kj and c = 2 A = V 2 , and obtain 
the Bethe-Ansatz equations for the multi-channel Ander- 
son impurity model. Removing the impurity {i.e. taking 
= 0) one recovers the usual equations of the NBA 
formalism. 66 



APPENDIX B: KERNELS AND IDENTITIES 



1 77 A 

K n (z) = (27ri)~ d z lne„ (z) = 



tt z 2 + {nAy 



plus the similar definitions 



K nm (z) = (2m) 1 d z In e nm (z) 
K ' nm 0) = (Zniy 1 d z lne' nm (z) 
K>; m (z) = (2 m y 1 d z lne'> m (z) . 

These last three definitions are the basis to define the 
following 'convolution' kernels: 

A nm ( z ) — S n , rn S (z) + K nm (zj — 

min{m,n} 

= K m+n+ i +a - 2 T (z) 

<7=±1 T=l 

min{m,n} 

B nm (z) = K' nm (z) = K m +n+i-2 T (z) 

r=l 

min{m,n} 

C nm (z) = 8 n . m 5 (z) + K" m (z) = ^ K m+n~2r (z) 



that are used extensively in the continuum formulation 
of the BAE. 



b. Fourier Space Formalism 

A great simplification in the algebraic manipulations is 
often achieved by working in terms of the Fourier trans- 
formed densities. Our convention will be as follows: 



In this appendix we collect a number of function and 
kernel definitions, and identities relating them, that are 
central to the writing, rewriting, and algebraic manipu- 
lation of the BAE and the equations of the TBA. Let us 
introduce the following notations aimed at lightening the 
writing of the BAE, 



e„ (2) = 



z — inA 
z + inA 

min{m.n} 



e 'nm ( Z ) — H e m+n+l-2T (z) 



T=l 

min{m,n} 



r=l 

Snm {z) = e nm+a ( Z ) ■ 



For the purpose of writing a continuum version of the 
BAE, we will need as well the derivatives of the loga- 
rithms of the above functions. We define the following 



p(w) = F{p{z)} (w) = / p(z) 



So that, for instance, the basic kernels adopt the simple 
form K n (w) = e~ nA ' w \. Working in Fourier space is 
easy to see that different convolution kernels are simply 
related: B n , m = KiC n ,m and B„ jm = GA n , m - Where 
one defines the 'basic recursion kernel', 



G(w) 



f 



K H + K 2 H 
G{z)=T- 1 |g(w)} (z) 



2 cosh Aw 
1 

4A cosh H 



that we call so because it enters many recursion rela- 
tions connecting the different convolution kernels (for 

instanced A n>m -GA n +i,m + {^ - S n ,i) GA n -i, m = $n, m 
plus many others involving also the convolution kernels 
B n ,m an( i G n , m ). All these relations are easy to prove in 
Fourier space. 
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General Recursion Kernels 



d. Rank Recursion Kernels 



The basic recursion kernel is ubiquitous in the TBA 
equations of all integrable models. For the case of the 
multi-channel Anderson model, we are also going to need 
the more general kernels: 

q(N,M) = ^ Nm = e -(N-M)A\w\ sinh (mAw) 
A MM sinh (M Aw) ' 

Of which the basic recursion kernel is a particular case, 

(2 2) 

G = G\ ' '. When m = N some of these kernels are 
singular (have non-zero asymptotics in Fourier space). 
We regularize them according to 

fl e M) (A) = G^ M ) (\)-S m . N S(X) 

where we have expressed them, in direct space, in terms 
of the variable A = ttz/2A. When carrying out numerical 
calculations, we will use the following explicit expressions 
for these kernels: 

Rfi(N,M) _ 

/ N + Tin a\\ 

er=±l x ' 

r=±l 

where F (z) = d z lnT (z) is the digamma function. 



We also need, for intermediate manipulations, kernel 
operators that act on the rank indices. The basic one 
is G rs = 5 r s - (5 r s +1 + V' 1 ) G. Using it we can extend 
the convolution kernels as A™ m = G rs A nm . A partic- 
ularly useful derived kernel is the one given by the in- 
verse: = [G _1 G sr ] , where the indices vary in the 
range 1, . . . , N — 1 or 1, . . . , M — 1 depending on the case 
(X = N,M). The explicit formula for this operator, in 
Fourier space, is 



K x - 



sinh [min (r, s) Aw] sinh [(X — max (r, s)) Aw] 
sinh (Aw) sinh (XAw) 



and two particularly useful cases are given by the general 
recursion kernels: 



(A) = (A) = 



1 . ir(M-r) 
ttM Sm M 



K 



nN-1 
N 



(A) = (A) 



1 



cos - -- M ■' +cosh^ 



cos 221 + cosh 2A 



(these will be used in the determination of the Schrieffer- 
Wolff limit). 
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Indices a and a transform according to the corresponding 
fundamental representation, the notation a indicates that 
the index transforms instead with the complex conjugate 
representation. 

The equality is in the sense of distributions [i.e. almost 
everywhere) . 

We remind the reader that the electron gas with lin- 
ear dispersion is highly degenerate. In other words, 
the two-electron free Hamiltonian hi,2 = —i(9i + 
d 2 ) has as a general solution the wave function 
exp(ifcixi + ik 2 x 2 ) [A aia2 0(xi - x 2 ) + (Si,2 A) ai a 2 6(x 2 - 
xi)], with Si ,2 arbitrary. As follows then from degenerate 
perturbation theory, one must choose already at zeroth- 
order the correct combination with respect to which to 
turn on a given perturbation. Our result corresponds to 
such a combination for the case when the perturbation is 
a multi-channel Anderson impurity. 

Notice that in order to achieve a formulation in terms of 
continuum distributions, the essential limit to be taken is 
N c — > 00. We remark that this is the same limit required 
in order to recover the field theory. 

See Ref. |4C|| for a proof that it is a good approximation 
to restrict ourselves to the s-channel or q-channel highest 
weight states in the case of SU (2); this proof can be ex- 
tended to the general SU (N) case. 

We introduce the notation S„.i that stands for (1 — 5„,i). 
Here fi is the total chemical potential acting on the elec- 
trons, defined as fx = fx c + [i s + fi q with fi s = i J2n=i 

and = jE!=i K™- 

The notation in the subindex means: an arbitrary n such 
that n < N. 

We remind the reader that, for the sake of simplicity, much 
of the discussion assumes /i = 0. 

Notice that T s 7^ Ts; these two scales only coincide, asymp- 
totically, in the quadrupolar moment limit (e s 2> e q ). 
Alternatively, the whole discussion can be rephrased in 
terms of e — fi. 

The different panels were nevertheless computed all inde- 
pendently as a check of the correctness of the program. 
Technically, in our case the SW transformation needs to 
be carried out carefully to the next order to see that the 
operators associated with the high-energy level mix with 
the low-energy ones and yield singular contributions when 
their associated degeneracy is the lower one. 
The differences stemming from the changed definition of 
the monodromy matrix: it will be the first index of the 
scattering matrices, rather than the second one, that shall 
be used as auxiliary space. 

For reasons of brevity we shall henceforth use either nota- 
tion: A a = A (a), u a = u (a), etc. 



